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Abstract. Gravastars have been recently proposed as potential alternatives to explain the 
astrophysical phenomenology traditionally associated to black holes, raising the question of 
whether the two objects can be distinguished at all. Leaving aside the debate about the 
processes that would lead to the formation of a gravastar and the astronomical evidence 
in their support, we here address two basic questions: Is a gravastar stable against generic 
perturbations? If stable, can an observer distinguish it from a black hole of the same mass? To 
answer these questions we construct a general class of gravastars and determine the conditions 
they must satisfy in order to exist as equilibrium solutions of the Einstein equations. For such 
models we peiform a systematic stabiUty analysis against axial-perturbations, computing the 
real and imaginary parts of the eigenfrequencies. Overall, we find that gravastars are stable to 
axial perturbations, but also that their quasi-normal modes differ from those of a black hole of 
the same mass and thus can be used to discern, beyond dispute, a gravastar from a black hole. 



PACS numbers: 04.40.Dg, 04.30.Nk, 04.25.Nx 
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1. Introduction 

The gravastar model was proposed by Mazur and Mottola f[\ and has attracted attention as 
a possible alternative to black holes. Within this model, in fact, a massive star in its late 
stages could end its life as a gravastar, namely a very compact object whose radius would be 
very close to the Schwarzschild radius (indeed arbitrarily close to it) without having an event 
horizon or a central singularity. For this to happen, a phase transition is expected to take place 
at or near the location where the event horizon would have been formed otherwise The 
interior of what would have been the black hole is replaced by a suitably chosen portion of 
de-Sitter spacetime with an equation of state (EOS) p = — p, surrounded by a thin shell of 
ultra-stiff matter with EOS p = +p, and which is then suitably matched to a Schwarzschild 
vacuum. In practice, the Mazur-Mottola gravastar model, hereafter the MM model HI > is a 
spherically symmetric and static five-layer solution of the Einstein equations, including two 
infinitesimally thin shells needed by the junction conditions of the metric O). 

For as much as it is ingenious, the gravastar model also challenges a building block in 
modern astronomy, namely the existence of astrophysical black holes. In spite of the many 
observational evidences in favour of black holes, it may indeed be fundamentally impossible 
to give an irrefutable observational proof for the existence of a black-hole horizon if only 
electromagnetic radiation is received |4j and thus tell apart a black hole from a gravastar, 
if the latter existed. Such a challenge has obviously attracted a lot of attention and several 
related models have been recently proposed in the attempt of providing answers to two basic 
questions: Once produced, is a gravastar stable to generic perturbations? If stable, can an 
external obverver distinguish a gravastar from a black hole of the same mass? 

In ref. O] it was argued that the MM solution is thermodynamically stable, but other 
stabiUty analysis are not so easy to perform because of the model's structure. Visser and 
Wiltshire Q have analyzed the radial stability of a simplified model with three layers and the 
stability was shown to hold for a number of configurations. This stability was generalized 
by Carter i6| for gravastar models with different exteriors. Other possibilities for the interior 
solution have also been considered: Bilic et al Q replaced the de-Sitter interior by a Born- 
Infeld phantom, while Lobo replaced the interior solution by one that is governed by the 
dark-energy EOS [8J, and Lobo and Arellano matched interior nonlinear electrodynamic 
geometries to the Schwarzschild exterior ||9l- 

To remove in part the complications produced by the infinitesimal shells in the MM 
model, Cattoen et al lITOl have found that fluid gravastars can be built if the fluid is confined 
to a given layer and has there anisotropic pressures. The latter essentially replace the surface 
tension which was introduced in the original model MM model by the matching of the metric 
in the infinitesimally thin shells. Although anisotropic and with rather arbitrary equations of 
state, these pressures have the appealing property of being continuous and thus of allowing 
one to build equilibrium models without the presence of infinitesimally thin shells and thus 
look more seriously into the issue of stability. Indeed, as a step to go beyond the construction 
of equilibrium models, DeBenedictis et al ifTDl made a first attempt to investigate the stability 
of gravastars through a qualitative analysis of axial perturbations. 

In order to provide more definitive answers to the questions mentioned above, we 
have constructed a general class of fluid gravastars with finite shell thickness and variable 
compactness. While this class is similar and has been inspired by the one proposed in ifTOlfTTI . 
it also differs in two important aspects. Firstly, while we also consider a fluid gravastar with 
anisotropic pressures, we try and reproduce the most salient features of the original MM model 
by creating an internal and an external region which reproduce a de-Sitter and a Schwarzschild 
spacetime at finite radii and not only asymptotically. Secondly, for these models we determine 
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precise bounds for the properties of the metric functions and the compactness of the gravastar 
that yield equihbrium solutions, thus restricting considerably the class of possible solutions. 

Using these models we have carried out a systematic investigation of the stability analysis 
of gravastars against axial -perturbations lfT2l[T3l . thus extending the results discussed in lfT4ll 
to our gravastar model and computing explicit eigenfrequencies. In this way we were able 
to conclude that: a gravastar is stable to axial perturbations and indeed it is possible to 
distinguish it from a black hole if gravitational radiation is produced. More specifically, 
we have found that for all the models considered the imaginary part of the eigenfrequencies 
is always negative, thus indicating stability against these perturbations. Furthermore, while 
it is always possible to build a gravastar with given compactness and thickness such that 
it will have the same oscillation frequency as that of a black hole with the same mass, 
the corresponding decaying time will be different. Our results thus provide a way of 
distinguishing observationally and beyond dispute a gravastar from a black hole. 

The paper is organized as follows: in section 2 we review the main features of the 
original MM model and obtain numerical solutions which enable us to impose bounds 
for the parameters of the MM solution. In section 3, instead, we present our model for 
a fluid gravastar with anisotropic pressures and discuss the bounds derived from general 
conditions imposed by the EOS and the properties of the metric functions. In Section 4 we 
outline the perturbation equations and the numerical methods employed in the study of axial 
perturbations of our model. In section 5 we discuss the results obtained for the QNMs and in 
section 6 we present our final conclusions. 



2. The MM gravastar model 

We first present a quick review on the main features of the original MM model m for they 
will be used in further developments in this paper and highlight that precise bounds exist for 
the existence of solutions even for the simplest gravastar model. 

In general, we consider a static and spherically symmetric line element 

ds2 = -/(r)d<2+ +r2df}2, (1) 
h{r) 

and the Einstein equations must be solved for a perfect fluid at rest, such that there are three 
different regions with the three different equations of state 

I. Interior: < r < ri , p = —p , 

II. Shell: ri<r <r2, p = +p, (2) 

III. Exterior: r2 < r , p ^ p = . 

In region I, p is a constant given by pv — SHq/Stt, and the metric is that of a de-Sitter 
spacetime, so 

f{r) = Ch{r) = C(l - H^r"^) , < r < n . (3) 

where C is an integration constant, whose value will be determined later [cf. eq. (O]. 

In region II, a new dimensionless variable w is defined as w = Snr'^p, in order to obtain 
the following set of equations 

dr dh 



1 — w — /i ' 



(4) 
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dh 
~h 

wf 



1 — w — h 
l + w-3h 



w 



= const. 



(5) 



(6) 



An analytical solution can be obtained in the case of a very thin shell, i.e., for ri r2. It is 
easy to see that (|4]l is satisfied if we take h = 1 ~ fh/r and Aifi{r) = 2dm{r) = 87rpr^dr^ = 
wdr. Therefore, in the thin-shell limit, (|4|i can be integrated to yield 



1 



m (1 + wY 



< 1 



(7) 



where e is an integration constant. From the continuity of the metric coefficients / and h at 
ri and r2, it can be shown that the integration constants e, C, M and i/o are given in terms of 
Ti, r2, wi and W2 by the relations 



ri 



, W2 

In — 



1 

W2 



1 

Wl 



C 



I +W2 
1+tVl 



1 



e(l + u'2)^ 

W2 



m{r2) ; 



(8) 



(9) 



(10) 



(11) 



Finally, in region III and because of Birkhoff's theorem, the spacetime is described by 
the Schwarzschild metric, 

2M 



f{r) = Kr) = 1 



r2 <r . 



(12) 



Clearly, region II is the most interesting one from a physical point of view since it is 
where a non-trivial model for the gravastar can be specified. In the case of a "thick shell", 
that is away from the thin-shell limit, the gravastar solution in region II can be obtained in 
the same way as for ordinary spherical stars, namely through the numerical solution of the 
Tolman-Oppenheimer-Volkoff (TOV) equations for p{r) = p{r) and m{r). We recall that in 
the MM model p{r) and p{r) are discontinuous in ri and r2 and that m{ri) must, of course, 
be less than ri/2 

By specifying ri, r2, and thus the shell thickness 6 = r2 — ri, as well as the initial 
conditions for the pressure p{ri) and mass m(ri), the numerical solution of the TOV 
equations provides m(r) and p(r) for ri < r < r2. A systematic analysis has revealed 
that a limit exists on the compactness /i = M / r2 of the gravastar and in particular that, for 
given ri, r2 and m(ri) very close to the limit ri/2, there exists a value for p{ri) which 
provides the largest M. In other words, it is not possible to achieve arbitrary values for M. 
On the other hand, if we fix r2 and decrease ri so as to increase the thickness of the shell, we 
must also decrease m(ri), and we have verified that the maximum value obtained for M also 
decreases. Overall, therefore, each shell thickness 6 also selects a maximum compactness of 
the gravastar p, above which no solution of the TOV equations can be found. 

All of this is shown in figure [T] which reports the space of parameters [5, p) where 
equilibrium solutions can be found. The solid curve, in particular, is computed numerically 
and distinguishes the region where equilibrium models can be found {i.e., the region below 
the curve), from the region where no solutions can be found {i.e., the region above the 
curve). Stated differently, for any given shell thickness 6, the solid curve marks the 
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maximum compactness {i.e., the largest value of /i) for which gravastar models can be 
built. Interestingly, therefore, it is not possible to build a very compact gravastar with a 
very thick matter shell. Rather, gravastars have either large compactness and thin shells, 
or small compactness and thicker shells, as shown in figure [T] This figure, and its inset in 
particular, also illustrates that gravastars can be built with arbitrarily large compactness, i.e., 
with /i ^ 1/2 and thus with the outer radius 7-2 being only infinitesimally larger than the 
corresponding Schwarzschild radius. It is exactly this property that makes gravastars hard to 
distinguish from a black hole if only electromagnetic radiation is available. 




Figure 1. Limits in tlie compactness and tliickness of a gravastar in the MM model. The curve 
shows the maximum compactness for a given thickness 5 of the shell, that is, in the light gray 
area below the curve we have the possible solutions, while in the area above the curve solutions 
are no longer possible. The dark gray area highlights the typical compactnesses for neutron 
stars. The inset shows a comparison between the numerical solution of the TOV equations 
with the analytical solution in the thin-shell limit. 



3. The gravastar model with anisotropic pressures 

3.1. The field equations 

We now discuss our suggestion for a fluid gravastar model in which we follow the spirit of 
the original MM model and consider a spacetime consisting of three different regions with 
the internal and external ones reproducing a de-Sitter and a Schwarzschild spacetime at finite 
radii, ri and r2, respectively. In addition, we follow the suggestion made in ifTOl and use 
a thick shell with a continuous profile of anisotropic pressures to avoid the introduction of 
infinitesimally thin shells of matter as demanded by the MM ||T|. We start therefore with a 
metric of the form 

d.s2 = -c'-t'-'dt^ + c^('-)dr2 + r'-An" , (13) 

and a fluid stress-energy tensor T^^^ ~ diag[~p, pr,pt,Pt], where pr andpt are the radial and 
tangential pressures, respectively. The Einstein field equations for this spacetime geometry 
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and matter distribution are 

e^M --- I +- = 87rp, (14) 

e"M - + - ) - - = S^Pr- , (15) 

It is now convenient to transform the above equations into a form where the hydrodynamical 
properties of the system are more evident and that reduces to the TOV equations for systems 
with isotropic pressure, i.e., 

e--^ = 1 \J- , (17) 

r 

, 2m(r) + ^irr^pr 



A' 

r 




1 




v' 
r 




1 




v" 


AV 




zy' - A' 


T 




4 


^ 2r 



r(r — 2m{r)) 
v' 

2+ . 



(18) 



Pr = -(P + PrjlT H , (19) 



where 

m(r) = / 47rr^pdr . (20) 
Jo 

Combining (fTSl ) and ( fT9t . we obtain the anisotropic TOV equation 

/ / , ,m(r)+47rrV , 2{pt - Pr) 

Pr = -{P+Pr)—, o f W + ' ^^^^ 

r[r — Zm[r)) r 
which is reminescent of the Newtonian hydrostatic-equilibrium equation and where the last 
term is obviously zero in the case of isotropic pressures, i.e., pt = Pr- 

3.2. Equation of state 

In order to adopt the model suggested in ifTOl . and still maintain the simple structure of the 
MM model, we make the following choices for our density function p{r) 

p(0) = p(ri) = po, P(r2) = 0, p'(ri) =p'(r2) =0. (22) 
A simple way to satisfy the above conditions is to consider a cubic dependence in r, 

{po , < r < ri region I. 

ar^ + br^ + cr + d , ri < r < r2 region II. , (23) 

, r2 < r region III. 

with the coefficients a, b, c, d given by 

2po 



(r2 - ri)3 



(24) 



, 3po ra +ri) 

b = -, , (25) 



(r2 - ri)3 ' 
(r2 - ri)3 



(26) 
(27) 
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To obtain the density in the terms of the total mass m(r2) = A/, we write 



(^2 - rif 
in 



(r2 - 3rir|)(r| - rj) rl{r2 - ri)^ 



(28) 

For the radial pressure Pr, we follow the suggestion made in refs. ifTSlfTTI and use an EOS of 
the type 

2" 

a - (1 



Prip) = — 
\P0 



(29) 



Clearly, the EOS (29i cannot be derived from basic principles and serves here essentially as 
a closure relation for the system of the equations. Yet, such an EOS can be constrained and 
the parameter a is determined by demanding that the maximum sound speed d^pr/dp^ = 
coincides with the speed of light to rule out a superluminal behaviour Because the sound 
speed 

2' 

a -2(1 



dpr 
dp 



Po 



.Po 



(30) 



has a maximum for p = p^^ aj (2(1 + a)), requiring Csifi) ~ 1 yields a ~ 2.2135. 
Despite being not particularly realistic, the EOS ( |29] | has the advantage of being simple and 
of possessing the needed physical limits since p^ —p for r ri and pr for r 7'2 
{cf. right panel of figure|2]for a typical gravastar model). 




Figure 2. Left panel: Behaviour of the functions p{r), Pr{i") and pt{r) for a representative 
gravastar model with M = 1, ri = 1.8 and r2 = 2.2. The function pt{r) is scaled by 0.1 for 
better visualization but provides the dominant contribution in eq. )2U . Right panel: Behaviour 
of the mass function m(r) and of the metric coefficient grr for the same solution shown in the 
left panel. 



Once the solutions for the energy density and the radial pressure are known, the 
tangential pressure pt can be computed through the anisotropic TOV equation (|2TI) as 

r 1 

Pt^Pr + -^P'r + -j{Vr^ P) 

The set of equations (I22l) -(|3T]) fully determines our gravastar model which will have a finite 
core described by the de-Sitter metric, a crust of matter and an exterior described by the 
Schwarzschild solution. 



m(r) + Arrr Pr 
7'(1 — 2m{r)/r) 



(31) 
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The behaviour of the energy density and of the pressures is shown in the left panel of 
figure|2]for a representative gravastar model with ]\I = 1, ri = 1.8 M and r2 = 2.2 M. The 
right panel on the same figure, instead, shows the mass function and the metric coefficient grr 
for the same representative model. 



3.3. Conditions on the metric frnctions 

While equations (l22l i-(l3Tli allow to construct gravastar models, additional constrains need to 
be imposed to guarantee that the metric functions gtt and grr (and their first derivatives) have 
the expected properties once our choice for the EOS Pr{p) and p{r) is made. In particular, it 
is not difficult to conclude that m(r) and m'{r) are continuous for the choice of p{r) made in 
equation ( |23] ), so that using equation (fTTI i grr and its first derivative are continuous throughout 
the spacetime. 

Similar considerations apply for the metric function gu, for which we note that in 
obtaining the solution for v by integrating equation ( fTSl l we must allow for a constant of 
integration vq so that 

^= / — 1 o f \\ dr + i/Q- (32) 

Jo r[r - 2m(r)) 

The integration constant can then be determined through the condition that, at the surface of 
the gravastar, gttir-i) = ^ (1 ^ 27\//r2), thus giving 

gu = -e'^t'-' = - (^1 - ^) ci^('-)-r(r.) ^ (33^ 

where 

m.r'^fl^^r. <34) 

Jo r{r - 2m{r)) 

In this way gtt and g'^^ are continuous across r2 and throughout the spacetime. 




-0.2 1 ' ' ' ' 1 -0.1 I ' ' ' ' 1 

0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 

r S 

Figure 3. Left panel: Behaviour of 1 — 2m{r) /r for M = 1 and different values of ri and 
r2 = 2.2. Right panel: Dependence on & of the minimum value of 1 — 2m(r) / r for gravastars 
with M = 1 and different values of r2 . Note that for sufficiently large models, the minimum 
is always positive. 

The requirement for the metric functions of being continuous with continuous first 
derivatives is, in general, not sufficient to guarantee the existence of acceptable equilibrium 
models. The three free parameters ri and r-z, in fact, cannot be chosen arbitrarily but in 
such a way that grr is always positive. In the left panel of figure[3]we show a typical example 
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of the behaviour of the function 1 — 2m{r)/r, where different lines refer to different values 
of the inner radius ri while the outer one r-2 is held fixed. Clearly, an incorrect choice of ri 
and r2 leads to negative values for g^,- and thus to the unphysical appearance of horizons. The 
right panel of figure[3]shows how the minimum value of the function 1 — 2m (r) / r depends on 
the gravastar's thickness 5 and that gravastars with large thicknesses can be built if the outer 
radius r2 is chosen to be sufficiently large, that is, if the compactness ^ is sufficiently small. 

As a more systematic characterization of this problem, we have constructed a large 
number of gravastar models in which we have varied both the compactness ji and the thickness 
6. In this way we can extend the diagram presented in figure [T| for the existence of MM 
models to our thick-shell gravastar model and define the region in the {5, i^i) parameter 
space where equilibrium solutions can be found. In particular, by definining the parameter 
e as the dimensionless distance of the gravastar's surface from a Schwarzschild horizon, 
i.e., e = r2/M — 2, we find that for e > ec « 0.3085, the function 1 — 2m{r)/r is always 
positive, so that the thickness 5 can be as large as r2 {i.e., ri can also be taken to be zero). 
On the other hand, for e < Ec there exists a critical thickness 5c = Sc{r2, M) above which 
equilibrium models cannot be found because of the appearance of horizons. For e ~ Ec, we 
have found that Sc = 2.30685 M. 




Figure 4. Limit on tlie compactness /i of tlie gravastar witli tlie thickness of tlie shell <5. This 
figure compares the results obtained for the MM model in the lower curve (see figure [D and 
the results obtained for our model in the upper curve. 

Figure |4] summarizes these results by showing the region of the parameter space {S, pi) 
where equilibrium models can be calculated. Equilibrium solutions for our model can be 
found in the region below the solid line, and 6/M = 2.30685 is the boundary between the 
two distinct behaviors discussed above. Indicated with a dashed line is the corresponding 
threshold for the MM model (cf. figure [Til. Overall, the data in the figure reveals that our 
thick-shell model obeys similar but less restrictive bounds than the MM model. 
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4. Perturbation equations and numerical methods 

In order to assess the stability of gravastars against axial perturbations we have followed the 
standard procedure to study non-radial oscillations of stars lfT2l [131 1161 [TtI . This approach 
differs from the one adopted in ref. lITTl . where the analysis of the axial perturbations of a 
gravastar was made following the standard procedure for Schwarzschild black holes (see lITSl 
for a review) in order to study only the asymptotic properties of the potential and determine 
whether it has compact support and is everywhere positive. 

We note that the pressure anisotropy in the shell does not play an important role in the 
study of axial tensor perturbations since these perturbations hardly excite fluid motions. As 
a result, the standard analysis for stars with isotropic pressures can be used here. This is no 
longer true for polar perturbations (for which the pressure anisotropy plays an important role) 
and a study on radial perturbations of stars with anisotropic pressures was recently presented 
in HDIlOl. 

As for Schwarschild black holes or spherical relativistic stars, the propagation of axial 
perturbations in the gravastar spacetime is governed by the wave equation in a scattering 
potential 

^-^I^^KO^, (35) 

where the "tortoise" coordinate is here defined as 

EE /" e'^-'^^/^dr, (36) 
Jo 

and the scattering potential is given by the expression 

Vi{r) EE ^ + l)r + 47rr3(p - p^) - 6to] . (37) 



The potential vanishes for r ^ oo and diverges as for r ^ 0, as a result of the centrifugal 
term proportional to + 1) jr^ . This is an important difference that the potential Viir) has 
with respect to its counterpart for a Schwarzschild black hole, where Vi(r) -^0 at the horizon 
T — 2 M. Note also that ey ^ \ for r ^ oo (the spacetime is the Schwarzschild one for 
r > Ti) and, from eqs. ^ and dH, e'^ ^ (1 - 2M/r2)e-^(''=) for r 0. 

Unlike for black-hole spacetimes, both the tortoise coordinate and the scattering potential 
are here not expressed through simple algebraic relations but need to be computed as solutions 
of ordinary differential equations. We have done so using a 4th order Runge-Kutta method to 
obtain <y and from equations ( l33l l and (l36l l, respectively. 

Introducing now the "light-cone" (null) variables u = t — and v = t + r^,, the wave 
equation (l35l l can be written in the more compact form as 

-i^^{u,v) = Veir)^{u,v), (38) 

and thus be integrated numerically. In particular, we have used a variation of the method used 
in II2TI and in many other works, and derived a strict second-order accurate discretization 
given by 

16 — A'^V 

= i^E + V'w) ,. , A,.; - V-s + OiA^) , (39) 



16 + 



s 



where the indices N, S, E and W refer to the grid-points defined as = (u + A, w + A), 
S = iu,v), E = {u,v + A), and W = {u + A,u). Equation ( l38l l was then numerically 
integrated in the (u, v) plane with the algorithm ( |39l l using a triangular grid limited by the 
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lines = r™*" and m = 0, where r™™ is a small value which we have set as r™™ = A/2. 
A schematic representation of the null grid is presented in fig. |5] where the black points 
represent the grid points where the solution is known, while the red ones are those where the 
solution is to be calculated. 



Figure 5. Diagram of the numerical grid and the domain of interest. The black points represent 
the grid points where the value of the field is known. The red points represent the grid points 
to be calculated. 



In a linear regime the eigenfrequencies of the gravastar are not sensitive to the choice of 
the initial data and thus we have modelled the initial perturbation with the simplest possible 
choice, i.e. a Gaussian pulse centred around Vc and with witdth a 



ip{u ~ uq, v) = cxp 



Vc) 
2cr2 



(40) 



where, to satisfy the regularity condition at the origin, we set 

^{rT"\t)=i;{u = v~vo,v)^0, Mt . (41) 

During the integration of eq. ( |38] ) we extract the values of the field ?/> along a line of 
constant r*, i.e., at the points {u — v — 2r,, v). We let the field evolve for large values of 
t, until the transient is over and possible contributions of the higher overtones have died off. 
Only then, when the signal consists of the fundamental mode (the slowest decaying mode) 
and of its overtones, we can compute the real oj^ and imaginary part lOj of its eigenfrequency 
with a least-squares fit of the function 

i;fu{t) = A„exp(co7i) cosiuj'lt + 4>) , (42) 

where A and are constant coefficients and n represents the order of the mode and this 
is not meant as an exponent. Once the properties of the fundamental mode are known, 
the eigenfrequencies of the overtones are obtained using the method proposed in 1221 . 
which consists of subtracting from the numerical solution the fitted function relative to the 
fundamental mode n = 0. This reveals the n = 1 mode and the procedure can be iterated 
for obtaining higher overtones, depending on their damping rate. Of course, very rapidly 
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decaying modes cannot be obtained by this procedure, but this approach is a very efficient 
and easy way to obtain the first overtones in a signal, which are the ones expected to be 
detected from any astrophysical sources. 

The numerical setup described above has been first tested by performing a perturbation 
analysis of Schwarzschild stars, i.e., uniform-density spherical stars, obtaining results that are 
in very good agreement with the values reported in the literature lfT3ll23l . 



5. Results 



Before discussing the results of the perturbative analysis, we note that the axial potential given 
in ( IJTI i shows features which are similar to those that are well-known in compact uniform- 
density stars ifTsl . Quite generically, in fact, the potential exhibits at least a minimum and 
a maximum, thus indicating the possibility of having modes that are trapped in the potential 
well. 

In figure |6] we present some typical examples of the scattering potential Vi{r) when 
i = 2 and which, in general, depends on both the compactness ^ (as for relativistic stars) and 
on the thickness of the shell 5. The left panel of figure |6] in particular, shows the potential 
of gravastar models having the same compactness ji = 0.46 but increasing thickness of the 
shell, while the right panel reports how the potential changes for gravastar models with the 
same shell-thickness 6 = 0.2 but increasing compactness. Clearly, in the exterior region Vt{r) 
is just the Schwarzschild potential for axial perturbations and the depth of the potential well in 
the interior region increases with both ^ and 6. Furthermore, given the chosen EOS (l29l l, it is 
also possible to select the parameters ri, r2, and M such that Vi{r) has two local minima very 
close to each other. This is shown by the curve with /^t = 0.35 in the right panel of figure|6] 




Figure 6. Left panel: The scattering potential Vi{r) witli I = 2of gravastar models with the 
same compactness ^ = 0.46 but increasing thickness of the shell. Right panel: The same as 
in the left panel but gravastar models with the same shell-thickness & = 0.2 but increasing 
compactness. 



Following the procedure outlined in the previous Section, we have integrated the 
perturbation equation dSSl l for a large variety of gravastar models differing in mass, 
compactness and thickness. Figure |2]reports some typical results and shows the time evolution 
of the perturbations for a representative gravastar with M = 1, and radii ri = 1.96, r2 = 2.26. 
The left panel, in particular, reports the total solution il}{t) as extracted at r* = 10 during 
the initial stages of the scattering off the potential and allows one to discern the "beating" 
caused by the superposition of overtones and other rapidly decaying components of the signal. 
Interestingly, this beating of the fundamental mode with overtones becomes more pronounced 
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Table 1. Some typical values obtained for the the first overtones of the gravastars QNM 
eigenfrequencies for axial perturbations with t = 2. The gravastars have M = 1,T2 = 2.26 
and <5 given in the Table. We also show for comparison the equivalent frequencies for a 
Schwarzschild black hole of the same mass and for a standard TOV star with uniform density 
with M = 1 and _R = 2.26M (note that for these stars one must have R/M > 2.25 QU). 



model 


n 


= 


n 


= 1 


n 


= 2 














(5 = 0.30 


0.3281 


2.481e-3 


0.4865 


6.264e-2 


0.6534 


1.590e-l 


S = 0.35 


0.2943 


7.081e-4 


0.4459 


3.202e-2 


0.5922 


1.093e-l 


& = 0.40 


0.2575 


1.543e-4 


0.4011 


1.227e-2 


0.5384 


5.814e-2 


Schwarzchild black hole 


0.3737 


8.896e-2 


0.3467 


2.739e-l 


0.3011 


4.783e-l 


Schwarzchild star 


0.1090 


1.239e-9 


0.1484 


3.950e-8 


0.1876 


5.470e-7 



with increasing /i and 5, underlining that the contribution of the overtones increases (i.e., the 
modes decay more slowly) with both /i and 5 {i.e., as the the depth of the potential well 
increases). A similar relation between the decay rate of the modes and the potential well is 
known for the "trapped modes" of compact uniform-density stars 123)1 . 




10 20 30 40 50 60 70 80 50 100 150 200 250 300 350 400 

1 t 



Figure 7. Left panel: Evolution of the axial perturbations tl){t) as extracted at r« = 10 
and during the initial stages of the scattering for a gravastar with M = 1, t\ = 1.96 and 
T2 = 2.26. Right panel: The same as in the left panel but when the contribution of the 
fundamental mode n = 0, is distinguished from that of the overtones n = 1,2 and 3. 

The right panel of figure |2l on the other hand, shows the evolution of the perturbations 
over a longer timescale and when the overtones have been "removed" from the general signal 
using the approach discussed in the previous Section. Note that, for increasing n, the period 
of the oscillations becomes smaller, and the decaying timescale shorter Using a least-squares 
fit of the numerical data to the ansatz (l42] | we have computed the real and imaginary parts of 
the QNM eigenfrequencies and the numerical results for some typical cases are presented in 
table [TJ where they are also compared to the corresponding values for a Schwarzschild black 
hole of the same mass and a Schwarzschild star of the same compactness. The tabulated 
values are all for an £ ^ 2 perturbation and refer to the fundamental mode and to the first 
two overtones since the third one does not usually have a sufficiently clear slope to allow for 
an accurate measurement. If necessary, an increase in spatial resolution would improve the 
decay of higher overtones and allow for their computation. 

Figure[8]is the most important of this paper and it reports the variation of the £ = 2, n = 
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QNM eigenfrequencies as a function of the parameters of the gravastar: /i and 6. The two 
panels report separately the behaviour of the real and imaginary parts and indicate that both 
the period of the oscillations = 2-k jio^ and the the damping time = 2t: \ increase 
with increasing 8 and [i. The two thick horizontal lines represent instead the corresponding 
frequencies for a Schwarzschild black hole of the same mass, i.e., Muj^ ~ 0.37367 and 
MlOj = —0.08896 1161 . Clearly, while it is always possible to select the thickness and 
compactness of the gravastar such that it will have the same oscillation frequency of a black 
hole with the same mass (c/, left panel), the corresponding decaying time will be different 
and about one order of magnitude larger (c/, right panel). Stated differently, a gravastar and 
a black hole with the same mass cannot have the same complex QNM eigenfrequency when 
subject to axial perturbations. This result, which we have here presented for the fundamental 
mode, can be shown to be true also for the first two overtones. 




Figure 8. Left panel: Behaviour of the real part o;^ of the £ = 2, n = QNM 
eigenfrequencies as functions of 5 for different values of fj,. Right panel: The same as in 
the left panel but for the imaginary part uij of the eigenfrequencies. In both panels, the thick 
hoiizontal lines represent the conesponding frequencies for a Schwarzschild black hole. 



6. Conclusions 

Although the gravastar model of Mazur and Mottola fT\ represents an ingenious solution 
of the Einstein equations in spherical symmetry, it has also challenged one of the most 
cherished foundations of modern astrophysics: i.e., the existence of astrophysical black holes. 
Gravastars, in fact, can be constructed to be arbitrarily compact, with an external surface 
which is only infinitesimally larger than the horizon of a black hole with the same mass. As a 
result, the electromagnetic emission from the surface of a gravastar will suffer of essentially 
the same gravitational redshift as that of a black hole, thus making it difficult, if possible at 
all, to distinguish the two when only electromagnetic radiation is available. 

Without entering the relevant debate about the physical processes that would lead to the 
formation of a gravastar or the astronomical evidence in support of their existence 1241 . we 
have here considered two more fundamental questions: Is a gravastar stable against generic 
perturbations? If so, can an external observer distinguish it from a black hole? The short 
answers to these questions are that: a gravastar is stable to axial perturbations and indeed it is 
possible to distinguish it from a black hole if gravitational radiation is produced. 

To reach the first of these conclusions we have constructed a general class of 
gravastar models that extends the one proposed by Mazur and Mottola by replacing the 
infinitesimal shell of matter with one having finite size 6 and variable compactness /x. These 
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equilibrium solutions of the Einstein equations have then been analyzed when subject to axial 
perturbations and the eigenfrequencies of the corresponding QNMs have been computed 
explicitely. For all of the cases considered, the imaginary part of the eigenfrequencies has 
always been found to be negative, thus indicating the stability of these objects with respect to 
this type of perturbations. 

To reach the second conclusion, instead, we have shown that the QNM spectra of a 
gravastar and that of a black hole of the same mass differ considerably. In particular, while it 
is always possible to select 5 and /i such that the gravastar has the same oscillation frequency 
as that of a black hole with the same mass, the corresponding decaying time will be different. 
As a result, the gravitational radiation produced by an oscillating gravastar can be used to 
distinguish it, beyond dispute, from a black hole of the same mass. 

We plan to extend our stability analysis also to polar perturbations and determine whether 
or not these intriguing objects possess modes of oscillation that do not have a counterpart in 
compact relativistic stars and may therefore hint to new solutions of the Einstein equations. 
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